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Pomerons and Lattices: a Progress Report 

UKQCD Collaboration - presented by C. Parrinello'^ 

'^Department of Mathematical Sciences, University of Liverpool, 
Liverpool L69 3BX, United Kingdom 

We report on some attempts to use lattice QCD to investigate topics in strong interaction phenomenology 
which are usually interpreted in terms of soft Pomeron exchange. 



1. INTRODUCTION 

It has been known for a long time that the rise 
of the total hadronic cross section at high energies 
could be described in the framework of Regge the- 
ory via Pomeron exchange 0. Renewed interest 
in Pomeron models has been recently triggered 
by the analysis of diffractive scattering at HERA. 
In particular, the discovery of events with large 
rapidity gaps has been interpreted as evidence 
for Pomeron exchange Recovering the phe- 
nomenology of Pomeron exchange from first prin- 
ciples QCD is a major theoretical problem. Here 
we describe some attempts to use lattice QCD to 
test models of diffraction and the Pomeron. We 
first describe an application of the method to the 
Landshoff-Nachtmann (LN) model Q and then 
we outline a more general approach to the lattice 
investigation of diffractive physics. 

2. THE LN POMERON 

In the LN model Pomeron exchange between 
quarks behaves like a C = -1-1 photon-exchange 
diagram, with amplitude 



iPl{u-i^,u){uY''u). 



(1) 



/3o represents the strength of the Pomeron cou- 
pling to quarks, and is related to the (non- 
perturbative) gluon propagator D{p) by 



/3o = 



1 



367r2 



(2) 



where g is the gluon-quark coupling. The model 
yields simple formulae for pp scattering, exclusive 
p production in deep inelastic scattering and the 



J 1'^— nucleon total cross section, which all con- 
tain integrals in momentum space oi g^D(jp) 

In order to make predictions one needs an ex- 
pression for g^D(p). By inserting in the model 
the gluon propagator as computed on the lat- 
tice (in the Landau gauge) we aimed to perform 
a strong consistency test from the point of view 
of QCD Q. For this purpose we generated two 
sets of quenched configurations and in addition 
we used data for the gluon propagator as eval- 
uated in [p|. The lattice parameters are listed 
in Table nT Full details of the results and the 
method for the determination of the gluon prop- 
agator can be found in ref. In the figures 
the results obtained from the data sets 1, 2 and 
3 are plotted with a dotted, solid and dashed line 
respectively. We inserted in the LN formulae an- 
alytical expressions corresponding to the best fits 
to the lattice gluon propagator on each data set. 
It turns out that our calculations are not signifi- 
cantly affected by discretisation effects, as all the 
relevant integrals in momentum space appeared 
to reach their asymptotic value at an integration 
scale well below the inverse lattice spacing. It was 
assumed that in the continuum limit the propa- 
gator is multiplicatively renormalisable, as it is 
in perturbation theory. Also, we neglected the 
running of the QCD coupling, i.e. we made the 
approximation g{p) = g- As the scale for the mo- 
menta in D(j)) was set from independent string 
tension measurements, we only had one free pa- 
rameter to fix in the expression g^D{p). This is a 
multiplicative factor corresponding to the prod- 
uct of a gluon wavefunction renormalisation con- 
stant times a numerical value for . We call this 
parameter gl^. It was determined by using (§ as 
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Table 1 

Parameters of the lattices used in the study of the LN model. 







Size 


Cfgs. 


a-' ( GeV) 


Data set 1 


6.0 


16" 


150 


1.9 


Data set 2 


6.2 




150 


2.7 


Data set 3 


6.0 


243 X 48 


500 


1.9 



a nonperturbative renormalisation condition, im- 
posing that (3q attains its phenomenological value 
of 2.0 GeV"^ ItII: 



f^o = ^J d'p bc'ff Aat(p)]' = 4 GeV- 



(3) 



We finally inserted g^ff^iat(p) in^ the formulae of 
the model. We adopted the analysis procedure of 
Halzen and collaborators B. 



2.1. Proton-proton elastic scattering 

The calculation of a^^^ and i.e. the energy- 
independent part of the total and elastic differ- 
ential cross section for proton-proton scattering, 
provides a benchmark for the LN model. Sin- 
gle Pomeron exchange is expected to dominate 
^ up to -t ~ 0.5 GeV^. The measured total 
and elastic differential cross sections are usually 
parametrised as follows: 



ftot 



da 
It 



X 0.168 , 

^] ^ (4) 

ml dt ^ ' 



For small t, the elastic differential cross section 
behaves like e^*, and the model is characterised 
by two parameters, a^^^ and B. 

We computed a'^^^^ and B on each lattice us- 
ing the lattice gluon propagator and the effec- 
tive coupling gcff- We obtained values of a^Q^ 
ranging from 18.12 mb to 19.85 mb and B from 
12.6 GeV"^ to 13.6 GeV"^, thus in very good 
agreement with each other, suggesting that both 
quantities are subject to only small finite vol- 
ume effects. They are also encouragingly close 
to the phenomenological values of a^^^ ~ 22.7 mb 
and B ~ 11 GeV"^ In Fig. |l| we show ISR 
data j|] for the differential elastic cross section at 
y/s ~ 53 GeV together with the lattice predic- 
tions, with the energy correction of Eq. H. 




Figure 1. Data for the pp elastic cross section at 
y/s — 53 GeV together with the lattice predic- 
tions, corrected for the energy dependence. 



2.2. J/i/;-nucleon scattering 

This process provides a further important test 
of the LN model. Two phenomenological features 
emerge from our calculations. Firstly, the quark- 
counting rule is closely satisfied for light mesons, 
as we get app/aTip in the range 1.5 — 1.8. Sec- 
ondly, the Pomeron couples more weakly to heav- 
ier mesons, as can be seen from Fig. |[ where we 
show the meson-nucleon cross section as a func- 
tion of the pole radius, together with a Regge fit 
to the energy-independent part of the n^p and 
K~p cross sections |l]. Sizeable discrepancies be- 
tween different lattices can be observed in this 
case, which may be related to finite volume ef- 
fects. In this connection, it is worth noting that 
for most quantities the results obtained from the 
largest lattice provide the best agreement with 
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Figure 2. The meson- nucleon total cross section 
as a function of the radius used in the meson form 
factor. Also shown are the radii and total cross 
sections for the tt, K and J/V'- 



the experimental values. In order to further in- 
vestigate this point we plan to repeat our calcu- 
lations for a wider range of lattice parameters. 

3. DIFFRACTION ON THE LATTICE 

Here we summarise a more general approach 
for lattice investigations of diffractive scattering. 
We start by assuming a factorisable Pomeron ex- 
change. To clarify this concept, consider diffrac- 
tive scattering at HERA in the one-photon ap- 
proximation: 



7*{q)+p{p) ~^p{p) + x{px)- 



(5) 



If the above process factorises, one can describe 
it in two steps : first the original proton emits a 
Pomeron, p{p) — > p(p) + lP{A), then the Pomeron 
collides with the 7* to produce the final hadronic 
state X, 7*(g) + JP(A) X{px). On the lat- 
tice we aim to study the first subprocess, calcu- 
lating the effective coupling of the Pomeron to 
the proton. One motivation for such a study is 
the investigation of the helicity structure of the 
Pomeron, as it has recently been argued that a 



nontrivial spin structure in the Pomeron-proton 
vertex may affect quantities measured at HERA 
1^ . Given some ansatz for a composite operator 
0]p{x) which creates the Pomeron from the QCD 
vacuum, the effective Pomeron-proton coupling 
takes the form of a QCD 3-point vertex function, 
which we can study in momentum space as a func- 
tion of the momentum of the proton, p, and the 
momentum of the Pomeron, A. We assume that 
Ojp only contains gluonic fields. As for the quan- 
tum numbers, since the Pomeron is a Regge tra- 
jectory all the information we have a priori is that 
Op should have a J++ structure. As we do not 
want to assign a value to J, it makes sense to al- 
low Op to have arbitrary Lorentz structure, i.e. 
to be the sum of a scalar, a vector, and tensors of 
arbitrary rank. 

Numerical work is in progress. We are cur- 
rently focusing on two-gluon operators for Op 
with encouraging preliminary results. 
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